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Abstract: 

The first calculation of the next-to-leading-order electroweak corrections to Z-boson + jet 
hadroproduction including leptonic Z-boson decays is presented, i.e. to the production of a 
charged lepton-anti-lepton final state in association with one hard jet at the LHC and the 
Tevatron. The Z-boson resonance is treated consistently using the complex-mass scheme, and all 
off-shell effects as well as the contributions of the intermediate photon are taken into account. 
The corresponding next-to-leading-order QCD corrections have also been recalculated. The full 
calculation is implemented in a flexible Monte Carlo code. Numerical results for cross sections 
and distributions of this Standard Model benchmark process are presented for the Tevatron and 
the LHC. 
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1 Introduction 



The Drell-Yan process is a cornerstone of electroweak (EW) physics at hadron colliders like the 
Tevatron and the LHC. The production of W and Z bosons (or off-shell photons) with subsequent 
leptonic decays has both clean signatures and large cross sections. It can potentially be used 
to measure the luminosity of the collider, to constrain the PDFs, or to calibrate the detector. 
On the theoretical side, it is a perfect testing ground for our understanding of hadron-collider 
physics. Combining experimental accuracy and theoretical predictivity allows for a number of 
precision measurements in spite of the hadron-collider environment. (See e.g. Refs. [HE] and 
references therein.) 

At hadron colliders, the EW gauge bosons are almost always produced together with addi- 
tional QCD radiation. In this work, we focus on the neutral-current Drell-Yan process in which 
the dilepton pair is produced in association with a hard, visible jet, i.e. 

PP/PP -> Z/7* + jet l + l~ +jet + X. (1.1) 

The production cross section of this process, which is widely dominated by resonant Z bosons, 
is large, and the final state of the hard-scattering process is completely reconstructable. When 
the transverse momentum of the jet is large, the dilepton pair will be boosted, and the process 
provides a source for high-energy lepton pairs with opposite charge (opposite-sign lepton pairs). 
The invariant-mass distribution of the leptons is dominated by the Z-boson resonance, and a 
good theoretical understanding of this resonance also for boosted Z bosons is a cornerstone of an 
efficient detector calibration. Providing high-energy lepton pairs and jet(s), Z+jet(s) production 
is not only a Standard Model (SM) candle process but also an important background for new- 
physics searches. Moreover, the process offers the possibility for precision tests of jet dynamics 
in QCD. 

The importance of Drell-Yan processes is also reflected in the effort to make the theoretical 
predictions as precise as possible. The differential cross section for W/Z production is known at 
next-to-next-to-leading order (NNLO) accuracy (i.e. at two loops) with respect to QCD correc- 
tions [3] and even up to N 3 LO in the soft-plus- virtual approximation [4j. The next-to-leading- 
order (NLO) QCD corrections have been matched with parton showers [5] and combined with a 
summation of soft-gluon radiation (see e.g. Ref. [6]), which is necessary to predict the transverse- 
momentum distribution of the EW bosons for small px- I n view of QCD only, the neutral- and 
charged-current Drell-Yan processes are very similar. However, concerning EW corrections, the 
production of W and Z bosons show considerable differences and have been addressed separately. 
The NLO EW corrections are known for the charged-current [THlOj as well as the neutral-current 
process [HHE] and the predictions have been refined in various ways, e.g. upon including multi- 
photon radiation [T0iri3HT6] , photon- induced processes [12 jll3pi6Ifl"8] , and EW corrections within 
the MSSM pHQS]. Also the interplay of QCD and EW effects has been investigated [19]. 

The cross sections for W/Z + ljet [20H22] and W/Z + 2jets [H] production at NLO QCD 
are known for a long time. Recently, NLO QCD results for W/Z + 3jets and even W + 4jets 
production (in leading-colour approximation) were presented [23J. W/Z + ljet production has 
also been matched with parton showers [24]. Moreover, approximate results are available for the 
NNLO QCD corrections to Z + jet/Z + 2jets production for observables with especially large 
X-factors |25j. In the EW sector, higher-order corrections to W + ljet production have been 
first analyzed in the on-shell approximation |26H28] and later extended to the full NLO EW 
corrections for the physical final state [29], i.e. a charged lepton, a neutrino, and a hard jet. 

As far as EW corrections to Z + ljet production are concerned, only the purely weak one-loop 
corrections in the SM have been investigated in the on-shell approximation [30,31], i.e. in this 
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calculation the Z boson is treated as a stable external particle and photonic corrections have been 
ignored. For Z bosons at large transverse momentum, requiring a large centre-of-mass energy, 
using on-shell Z bosons is a good approximation since the EW corrections are dominated by 
large universal Sudakov logarithms [32]. In Ref. [30] the leading corrections up to the next-to- 
leading logarithms at the one- and two-loop level have been calculated. Later the full NLO weak 
corrections have been added [31]. However, the on-shell calculation is limited to a particular 
kinematic regime where neither off-shell effects nor the event definition for the physical final 
state play a role. Finally, photonic corrections have to be taken into account for a precision at 
the level of a couple of percent. 

In this work, we present a calculation of the NLO (i.e. one-loop) EW corrections for the 
physical final state, i.e. for the process pp/pp -4 / + /~+jet + X. Following our earlier work [29] on 
W+l jet production, the Z-boson resonance is described in the complex- mass scheme [33,34J. All 
off-shell effects due to the finite width of the Z boson, the contributions of and the interference 
with an intermediate photon, and photonic corrections are included. Our results have been 
implemented in a fully flexible Monte Carlo code which is able to calculate binned distributions 
for all physically relevant Z + 1 jet observables. In real emission events with photons inside a jet, 
we distinguish Z + jet and Z + photon production by a cut on the photon energy fraction inside 
the jet employing a measured quark-to-photon fragmentation function [35] . 

Our calculation is completely generic and thus not limited to specific observables or kinematic 
regimes. The interplay between the (potentially resonant) Z boson and the off-shell photon is 
included without approximations at NLO and can also be studied for observables for which the 
exchanged Z boson is far off shell and the photon contribution is potentially sizeable. However, 
we rely on the presence of a hard jet with sizable transverse momentum in the final state. For 
final-state jets which become soft or collinear to the beam pipe, the calculation breaks down. 
Nevertheless, the calculation of the EW corrections for Z production in association with a hard 
jet is also a step towards the mixed NNLO EW and QCD corrections to inclusive Z-boson 
production. 

We have also recalculated the NLO QCD corrections at 0{a 2 a1) in a fully flexible way, 
supporting a phase-space dependent choice for the factorization and renormalization scales, as 
will be discussed in some detail in Section [3l 

This paper is organized as follows. In Section [21 we describe our calculation in detail and 
discuss all the theoretical concepts and tools which have been used. In Section [3J we specify 
the numerical input as well as the details of our event selection. Numerical results are given for 
Z + jet production both at the LHC and at the Tevatron. We present inclusive cross sections for 
specified sets of cuts as well as distributions for the relevant observables. Our conclusions are 
given in Section |H 

2 Details of the calculation 

2.1 General setup 

At hadron colliders, the production of a charged-lepton pair via Z-boson or photon exchange in 
association with one hard jet is governed at leading order (LO) by quark-antiquark fusion, where 
the initial-state quarks radiate a gluon, and the corresponding crossed channels with a gluon in 
the initial state. To be specific, the relevant partonic processes are 




(2.1) 
(2.2) 
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Figure 1: Feynman diagrams for the LO process (|2.1]h 



qi g -> Z/ 7 *gi -> /+/" $ , (2.3) 

where % denotes any light quark, i.e. = u, d, c, s, b. The corresponding tree-level Feynman dia- 
grams for process (|2.1[) are shown in Figure [TJ The intermediate Z-boson resonance is described 
by a complex Z-boson mass /iz via the replacement 

M| M l = Ml - iM z r z (2.4) 

in the Z propagator as dictated by the complex-mass scheme (see below). Hence, all our results 
correspond to a fixed-width description of the Breit-Wigner resonance. Moreover, all related 
quantities, in particular the weak mixing angle, have to be formulated in terms of the complex 
mass parameters. The final-state leptons are treated as massless unless their small masses are 
used to regularize a collinear divergence. 

The tree-level amplitudes do not depend on the quark generation and only differ for up- 
and down-type quarks due to the different quantum numbers in the Zqq vertex. Hence, the 
summation over the quark flavours is straightforward for each of the three process types shown 
in (|2 . 1 [) — f|2 . 3 [) when folding the squared tree-level amplitudes with the corresponding PDFs. The 
five quark flavours (including the bottom quark), which appear as external particles, are treated 
as massless throughout the calculation, except if small masses are used to regularize a collinear 
divergence. At tree level the bottom-quark-induced processes do not show any peculiarities. 
Only in the evaluation of the EW virtual corrections they have to be treated with special care 
(see Section [272]) . 

In this work, we describe Z + jet production at NLO accuracy w.r.t. EW contributions, i.e. 
at the order 0(a 3 a s ). Hence, we also include the tree- level processes with a photon in the initial 
state, 

Ql 7 -)• Z/ 7 * qi -> l+r Ql , (2.5) 
q t 7 -> Z/ 7 *% -> l+r q % , (2.6) 

which contribute at the order 0(a 3 ) and may thus lead to relevant corrections at the expected 
accuracy level of a few percent. The tree- level Feynman diagrams for process (|2.5p are shown 
in Figure [2j The photon content of the proton has been quantified in the MRSTQED2004 
PDFs [36J. Since the photon also couples to the charged leptons in the final state, the amplitude 
is more involved than its QCD counterpart. As in our earlier work on W + jet production, in 
this work, we do not consider the crossed processes corresponding to l + l~ + photon production, 
which would lead to tiny corrections. The non-trivial definition and separation of the + jet 
and l + l~ + photon final states when additional photons are present due to bremsstrahlung are 
discussed in detail in Section 12.31 

To complete the description of the general setup of our calculation, we also repeat a few 
points which do not differ from our earlier calculation for W + jet production. To define the 
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Figure 2: Feynman diagrams for the photon-induced process (|2.5p . 



electromagnetic coupling constant a, we use the scheme, i.e. we derive a from the Fermi 
constant according to 



In this scheme, the weak corrections to muon decay Ar are included in the charge renormal- 
ization constant (see e.g. Ref. [8j). As a consequence, the EW corrections are independent of 
logarithms of the light-quark masses. Moreover, this definition effectively resums the contribu- 
tions associated with the running of a from zero to the weak scale and absorbs some leading 
universal corrections oc G^m^ from the p parameter into the LO amplitude. 

For corrections due to collinear final-state radiation it would be more appropriate to use a(0) 
defined in the Thomson limit to describe the corresponding coupling. On the other hand, using 
«G M everywhere is best suited to describe the large corrections due to Sudakov logarithms in 
the high-energy regime. Thus, the optimal choice cannot be achieved in one particular input 
scheme, and necessarily the calculation requires more refinements beyond NLO. In particular, 
among other things, higher-order effects from multi-photon emission should also be included at 
this level of precision which is beyond the scope of this work. We find that the difference of the 
two schemes in an NLO calculation only amounts to about 3% of the EW corrections. 

We employ the traditional Feynman-diagrammatic approach to calculate all relevant ampli- 
tudes in the 't Hooft-Feynman gauge. For a numerical evaluation at the amplitude level we use 
the Weyl-van-der-Waerden spinor formalism. To ensure the correctness of the presented results 
we have performed two independent calculations which are in mutual agreement. 

One calculation starts from diagrammatic expressions for the one-loop corrections generated 
by FeynArts 1.0 [37] • The algebraic evaluation of the loop amplitudes is performed with 
an in-house program written in Mathematica, and the results are automatically transferred to 
Fortran. The Born and bremsstrahlung amplitudes are calculated and optimized by hand and 
directly included into a Fortran program for numerical evaluation. A specific parametrization of 
phase space is used for an adaptive Monte Carlo integration employing the Vegas [38] algorithm. 

The second calculation is based on FeynArts 3.2 [39] and FormCalc version 3.1 [3D] . The 
translation of the amplitudes into the Weyl-van-der-Waerden formalism as presented in Ref. [H] 
is performed with the program Pole [42]. Pole also provides an interface to the multi-channel 
phase-space integrator Lusifer [33] which has been extended to use Vegas in order to optimize 
each phase-space mapping. Madgraph [34] has been very useful for internal checks of the 
real-emission amplitudes. 

2.2 Virtual corrections 

We calculate the virtual one- loop QCD and EW corrections for the partonic processes (I2.1|) - (l2.3p . 
but do not include the NLO QCD corrections to the photon-induced processes which are formally 




(2.7) 
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part of the corrections up to 0(a 3 a s ). The cross section of the photon- induced processes turns 
out to be numerically small at LO, and the corresponding corrections in the case of W + jet 
production turned out to be completely negligible. Hence, we do not expect phenomenologically 
relevant corrections for Z + jet either. We do not include the (loop-induced) contributions of the 
partonic process gg — > Z g. Based on Ref. [22J its contribution at the LHC and the Tevatron can 
be estimated to be below one percent. 

The calculated virtual QCD corrections are straightforward to implement and consist of up 
to box (4-point) diagrams only. The NLO EW corrections are more involved and additionally 
include pentagon (5-point) diagrams. There are 0(200) diagrams per partonic channel, including 
9 pentagons and 32 boxes. The general structure of the contributions is completely equivalent to 
our earlier work on W + jet production. However, Z + jet is computationally more demanding, 
not only because there are more diagrams per partonic channel but also because there are more 
helicity combinations contributing to the cross section which are summed at each phase-space 
point. The generic structure of the contributing diagrams is indicated in Figure El and the 
pentagon diagrams are explicitly given in Figure [U The different channels are related by crossing 
symmetry. 

Concerning the EW corrections, the partonic processes with one or two (anti-) bottom quarks 
in the initial state play a special role: There are box and pentagon diagrams with two W bosons 
and a heavy top quark in the loop. For all other partonic processes, only massless quarks 
propagate in the loops because we neglect CKM mixing, i.e. we set the CKM matrix to unity in 
our calculation. The mass of the top quark shifts the total EW correction roughly by a per mille 
for the most inclusive cross section discussed in Section El 

The potentially resonant Z bosons require a proper inclusion of the finite gauge-boson width 
in the propagators. We use the complex-mass scheme, which was introduced in Ref. [33] for 
LO calculations and generalized to the one-loop level in Ref. [34] • In this approach the W- 
and Z-boson masses are consistently treated as complex quantities, defined as the locations of 
the propagator poles in the complex plane. This leads to complex couplings and, in particular, 
a complex weak mixing angle. The scheme fully respects all relations that follow from gauge 
invariance. A brief description of the complex-mass scheme can also be found in Ref. [45] . 

The amplitudes can be expressed in terms of "standard matrix elements" [46], which com- 
prise all polarization-dependent quantities (spinor chains, polarization vectors), and invariant 
coefficients, which contain the tensor integrals. The tensor integrals are recursively reduced to 
master integrals at the numerical level. The standard scalar integrals are evaluated using two 
in-house Fortran libraries which are based on the results of Ref. |47] for 4-point functions with 
complex internal masses and on the methods and results of Ref. [38] . Results for different regu- 
larization schemes are translated into each other with the method of Ref. [49] . Tensor and scalar 
5-point functions are directly expressed in terms of 4-point integrals [50-52J , while tensor 4-point 
and 3-point integrals are recursively reduced to scalar integrals with the Passarino-Veltman al- 
gorithm [53]. Although we already find sufficient numerical stability with this procedure, we 
additionally apply the dedicated expansion methods of Ref. [52] in exceptional phase-space re- 
gions where small Gram determinants appear. 

UV divergences are regularized dimensionally. For the infrared (IR), i.e. soft or collinear, 
divergences we either use pure dimensional regularization with massless gluons, photons, and 
fermions (except for the top quark), or alternatively pure mass regularization with infinitesimal 
photon, gluon, and small fermion masses, which are only kept in the mass-singular logarithms. 
When using dimensional regularization, the rational terms of IR origin are treated as described 
in Appendix A of Ref. [54]; the ones of UV origin are always automatically included in the 
amplitude and tensor integral reduction. 
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Self-energy insertions: 




Triangle insertions: 





Box and pentagon insertions: 
I 




Figure 3: Contributions of different one-particle irreducible vertex functions (indicated 
as blobs) to the LO process ()2.2|) ; there are contributions from self-energies, triangles, 
boxes, and pentagon graphs. 




Figure 4: Virtual pentagon contributions to the process (|2.2[) , Note that for external 
bottom quarks the exchange of two W bosons leads to diagrams with massive top-quark 
lines (q' = top) in the loop. 
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We use an on-shell renormalization prescription for the EW part of the SM as detailed in 
Ref. |34| for the complex-mass scheme. Employing the scheme for the definition of the fine- 
structure constant a, we include Ar in the charge renormalization constant as mentioned above. 
The strong coupling constant is renormalized in the MS scheme with five active flavours. Hence, 
bottom quarks are included everywhere in the calculation as a massless quark flavour. 

2.3 Real corrections 

The evaluation of the real corrections has to be done with particular care, both for theoretical 
consistency as well as to match the experimental observables as closely as possible. Let us first 
focus on the EW real corrections to the partonic processes fj2.1|) — f)2.3[> . The emission of an 
additional photon leads to the processes 



The Feynman diagrams contributing to the process ()2.9[) are shown in Figure [5j Due to the 
emission of soft photons the real corrections include soft singularities which are cancelled by 
corresponding contributions in the virtual corrections independently of the details of the event 
selection or recombination procedure. If the photon and the charged leptons/quarks are recom- 
bined into a pseudo-particle (mimicking the start of hadronic or electromagnetic showers) to form 
IR-safe observables, all the remaining singularities arising from collinear photon emission in the 
final state also cancel against the corresponding singularities in the virtual corrections. This re- 
quires that all the selection cuts for a given observable are blind to the distribution of momenta 
in collinear lepton-photon configurations. The left-over collinear singularities due to collinear 
photon emission off the initial-state quarks are absorbed by a redefinition of the PDFs. Techni- 
cally, we use the dipole subtraction formalism as specified for photon emission in Refs. [55P56] 
to isolate all the divergences and ensure the numerical cancellation. 

Note that only the MRSTQED2004 PDFs properly account for all QED effects at NLO. 
However, this PDF set is outdated by now and does not include many modern PDF developments. 
Hence, we only employ the MRSTQED2004 set to estimate the photon content of the proton 
and use more modern sets for all the partonic channels without photons in the initial state. 
In a strict sense, with the choice of a modern PDF set the calculation is not fully consistent 
(like a NLO QCD calculation employing a LO PDF set). However, the numerical effect from a 
proper inclusion of the NLO EW corrections in the PDF determination is expected to be small, 
since this was also the case for the MRSTQED2004 fit when compared to the corresponding fit 
neglecting QED effects. The MRSTQED2004 PDF was defined in the DIS scheme so we stick to 
this scheme in our calculation. 

For muons in the final state it is experimentally possible to separate collinear photons from 
the lepton, i.e. to observe so-called "bare" muons. On the theoretical side, this corresponds to the 
fact that the lepton mass cuts off the collinear divergence in a physically meaningful way. Hence, 
the corresponding collinear singularities show up as logarithms of the small lepton (muon) mass 
since the KLN theorem [57] does not apply to non-collinear-safe observables. Our treatment and 
the analytical extraction of these collinear mass logarithms using the algorithm of Ref. [56] has 
been explained in detail in Ref. [29]. For Z + jet production, there are of course two charged 
leptons in the final state which can emit collinear photons. Apart from this straightforward 
generalization, which simply leads to more dipole subtraction terms, the same formalism applies. 



Qi Qi -> / + /~g7, 
qig -f l + l~ qi 7 , 
q~i g ->■ l + r q~i 7 . 



(2.8) 
(2.9) 
(2.10) 
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As in all processes with jets in the final state, the inclusion of EW corrections to Z + jet pro- 
duction asks for a precise event definition in order to distinguish Z + jet from Z + 7 production. 
We follow the strategy used for W+jet production which is detailed in Refs. [29,58]: We exclude 
jets which primarily consist of a hard photon (see Section [3]) and capture the non-perturbative 
physics in the collinear quark-photon splittings [59] by means of the measured fragmentation 
function [35] , Note that our calculation also provides an important step towards the inclusive 
prediction of dilepton pairs at large transverse momentum up to C(a 3 a s ), i.e. dilepton pairs 
recoiling against a jet or a photon, which would be free of the need to distinguish a photon from 
a jet [28]. The only missing ingredient is the combination of our results with the NLO QCD cor- 
rection for the dilepton+photon final state [60l[61]. However, since dilepton+photon production 
is suppressed with respect to the dilepton+jet final state, the inclusive results will at most differ 
at the level of one percent from the results presented in this work. 

Concerning the real corrections in NLO QCD, there are no particular complications when all 
the contributing partonic channels are identified. An additional gluon leads to the processes 

gift -►rT gg, (2.11) 

qig^l+l-qig, (2.12) 
Qig^l+l-qig, (2.13) 
gg ^l+l-q iqi . (2.14) 

Furthermore, the gluon present at LO may split into two quarks, inducing the processes 



qi qi ^ t 


+ l qi qi , 




(2.15) 


q%q% ->• 1 


+r qj qj 


(ft ^qj), 


(2.16) 


qi qj ->• i 


+r qi qj 


(qi ^qj), 


(2.17) 


q~i q~i -> I 


+ l~ q~i q~i , 




(2.18) 


qi q~j -> > 


+ r qi qj 


(qi ^ 1 


(2.19) 


qi q% ->• / 


+ l~ qi qi , 




(2.20) 


qi qj ->• i 


+ r qi qj 


(ft + qj) 1 


(2.21) 
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where qi denotes again quarks of the five light flavours. Note that different Feynman diagrams 
contribute for i = j and i 7^ j, so we list the corresponding partonic processes separately. Tak- 
ing this difference into account, as well as the correct quantum numbers for up- and down- type 
quarks, the remaining sums over flavour can again efficiently be performed when convoluting 
the squared matrix elements with PDFs. As in the EW case, we use the dipole subtraction 
method [62] to extract the IR singularities analytically from the numerical phase-space inte- 
gration. Absorbing all the collinear singularities due to initial-state splittings into the relevant 
PDFs, the remaining collinear and soft divergences cancel all the divergences of the one-loop QCD 
corrections for processes (|2.ip - (|2.3p . As explained already in Section [231 we do not include the 
NLO QCD corrections to the photon-induced processes. 

For the six-fermion processes (|2.15p . (|2.18|) . and (|2.20p with identical fermions, diagrams 
with gluon exchange can interfere with purely EW diagrams at C(a 3 a s ). The result is non- 
singular in the collinear limits due to the restrictions from colour flow, but in contrast to the 
other subprocesses with different quark flavours it does not vanish. However, the effect of this 
non-trivial interference contribution is expected to be phenomenologically negligible, as shown 
for the similar process W + jet production in Ref. [29] . Therefore, we neglect these contributions 
in our numerical analysis. 

3 Numerical results 

3.1 Input parameters and setup 

The relevant SM input parameters are 

Gf, = 1.16637 x Kr 5 GeV~ 2 , a s (M z ) = 0.1202, 
M$ s = 80.398 GeV, = 2.141 GeV, 

Mg s = 91.1876 GeV, T% s = 2.4952 GeV, M H = 120 GeV, 

m e = 0.510998910 MeV, m M = 105.658367 MeV, m t =172.6GeV, (3.1) 

which essentially follow Ref. [63]. As stated before, the CKM matrix only appears in loops and 
is set to unity, because its effect is negligible there. 

Using the complex-mass scheme [34], we employ a fixed width in the resonant W- and Z-boson 
propagators in contrast to the approach used at LEP and Tevatron to fit the W and Z resonances, 
where running widths are taken. Therefore, we have to convert the "on-shell" (OS) values of 
My S and Ty S (V = W, Z), resulting from LEP and Tevatron, to the "pole values" denoted by 
My and Ty. The relation between the two sets of values is given by [64] 

m v = M° s /Vi + (r£ s /M° s ) 2 , T V = r£ s /Vi + (r? s /M£ s ) 2 , (3 . 2 ) 

leading to 

M w = 80.370 . . . GeV, T w = 2.1402 . . . GeV, 

M z = 91.153 . . . GeV, T z = 2.4943 . . . GeV. (3.3) 

We make use of these mass and width parameters in the numerics discussed below, although the 
difference between using My or My S would be hardly visible. 

As explained in Section f2.il we adopt the G^ scheme, where the electromagnetic coupling 
a is set to olg^- In this scheme the electric-charge renormalization constant does not contain 
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logarithms of the light-fermion masses, in contrast to the a(0) scheme, so that the results become 
practically independent of the light-quark masses. 

We use the central MSTW2008NLO PDF set [65J in its LHAPDF implementation [66] for 
our numerical results. This implies the value of a s (Mz) stated in (|3.ip . We use the a s running as 
provided by the LHAPDF collaboration. Only the photon-induced processes are evaluated with 
the MRSTQED2004 set of PDFs [36] as discussed in Section 12.31 Here, we use the corresponding 
value of a s (M z ) = 0.1190. 

The QCD and QED factorization scales as well as the renormalization scale are always iden- 
tified. For low-px jets, the scale of the process is given by the invariant mass of the leptons 
which in turn peaks around M% for resonant Z-boson production. Hence, one natural choice is 
the Z-boson mass, i.e. /xr = = M%. For high-pr jets, well beyond the Z-mass scale, however, 
the relevant scale is certainly larger, and the QCD emission from the initial state is best mod- 
elled by the pt of the jet itself (see e.g. Ref. [67]). To interpolate between the two regimes, we 
alternatively use 



as a phase-space dependent scale, where p^ is given by the pt of the summed four-momenta of 
all partons, i.e. quarks and/or gluons, in the final state. At LO, p^ d is simply the pr of the one 
final-state jet. We present numerical results for both scale choices. 

3.2 Phase-space cuts and event selection 

In order to define IR-safe observables for the process pp/pp — > Z/7* + jet — > l + l~ + jet + X we 
recombine final-state partons and photons to pseudo-particles and impose a set of phase-space 
cuts as detailed in the following subsections. 

3.2.1 Recombination 

To define the recombination procedure and the separation cuts, we use the variables Rij = 
\J {m ~ Uj) 2 + (pfji where yi denotes the rapidity y = \ \a\{E + pi J )/(E — p^)] of particle i and 4>ij 
is the azimuthal angle in the transverse plane between the particles i and j. In the definition of 
the rapidity, E denotes the particle's energy and pi, its three-momentum along the beam axis. 
The recombination procedure, where we simply add four-momenta to form a pseudo-particle, 
works as follows: 

1. For observables with bare muons we do not recombine photons and leptons. For inclusive 
observables, a photon and a lepton are recombined for i? 7 / < 0.1. If both charged leptons 
in the final state are close to the photon we recombine it with the lepton leading to the 
smallest i? 7 /. 

2. A photon and a parton a (quark or gluon) are recombined for R^ a < 0.5. In this case, we 
use the energy fraction of the photon inside the jet, 27 = E^/(E^ + E a ), to distinguish 
between Z + jet and Z + 7 production. If 27 > 0.7, the event is regarded as a part of Z + 7 
production and rejected because it lacks any other hard jet at NLO. This event definition is 
not collinear safe and requires the use of quark-to-photon fragmentation functions to include 
the non-perturbative part of the quark-photon splitting as mentioned in Section 12.31 Our 
results are not very sensitive to the specific choice of the cut on 27. 




(3.4) 
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3. Two partons a, b are recombined for R a b < 0.5. For our simple final-state configurations, 
this procedure is equivalent to the Tevatron Run II /cp-algorithm [68J for jet reconstruction 
with resolution parameter D = 0.5. 

Technically, we perform a possible photon-lepton recombination before the photon-parton 
recombination. This procedure is IR safe because the triple-soft/collinear situation that a photon 
should have been first recombined with a parton, but was erroneously first recombined with a 
lepton, is excluded by our basic cuts. 

3.2.2 Basic cuts 

After applying the recombination procedure of the previous section we define Z + jet events by 
the following basic cuts: 

1. A partonic object (after a possible recombination) is called a jet if its transverse momentum 
Pt is larger than p™J et = 25 GeV. Events are required to include at least one jet. 

2. We demand two charged leptons with transverse momenta > 25 GeV. 

3. For the dilepton invariant mass we require Mu > 50 GeV to cut events with nearly on-shell 
photons splitting into a collinear lepton pair. 

4. The events have to be central, i.e. the leptons and at least one jet have to be produced in 
the rapidity range \y\ < y max = 2.5. 

5. The leptons have to be isolated, i.e. the event is discarded if the distance between one of 
the leptons and a jet Rij e t is smaller than 0.5. The lepton-jet separation is also required 
for jets with \y\ > y max . 

Note in addition that it is important to exclude low-px partons from the lepton-jet sep- 
aration procedure (guaranteed by step 1.), since otherwise observables would not be IR 
safe. 

While the EW corrections differ for final-state electrons and muons without photon recom- 
bination, the corrections become universal in the presence of photon recombination, since the 
lepton- mass logarithms cancel in this case, in accordance with the KLN theorem. Numerical 
results are presented both for photon recombination and for bare muons. 

For certain observables, we apply a jet veto against a second hard jet. To be specific, we veto 
any sub-leading jet with > Pr,ji/2, where Pt,j'i denotes the pt of the "leading" jet, i.e. the 
one with maximal p^. 

3.3 Results on cross sections and distributions 

We consider the production of a lepton pair in association with a jet at the Tevatron, i.e. for a 
pp initial state with a centre-of-mass (CM) energy of yfs = 1.96 TeV, and at the LHC, i.e. for a 
pp initial state. For the latter, we show results for y/s = 7 TeV, corresponding to the available 
energy in the years 2010 to 2012, as well as y/s = 14 TeV, the ultimate energy reach of the LHC. 

We present the LO cross section c>"o and various types of corrections 5, defined relative to 
the LO cross section by a = <7o x (1 + 5). Concerning the EW corrections, we distinguish 
the cross section for bare muons and Opyw for which a lepton— photon recombination is 

employed as defined above. Accordingly, the corresponding corrections are labelled d^yj and 
^EW' respectively. An additional label specifies which renormalization and factorization scale 
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is used. Either we use the fixed scale (/i = Mz) or we determine the scale on an event-by- 
event basis by the kinematical configuration of the final state (var), as specified in (|3.4p . For 
the EW corrections the difference is small, since the LO and the NLO results depend on the 
renormalization scale for a s and the QCD factorization scale in the same way. However, for 
the QCD part a sensible scale choice can be crucial for the stability of the perturbative series. 
Accordingly, the QCD corrections are labelled ^qcrf 2 f° r a fi xe d scale choice and 5qq^ for the 
scale choice defined in (|3.4p . 

As already observed for W+jet production in Ref. [29J and shown below, the QCD corrections 
become larger and larger with increasing of the leading jet. The increase in the cross section 
results from a new kinematical configuration which is available for the Z + 2 jets final state. The 
large of the leading jet is not balanced by the leptons, as required at LO, but by the second jet. 
Hence, we encounter the production of 2 jets where one of the quark lines radiates a relatively 
soft Z boson or off-shell photon. This part of the cross section, which does not really correspond 
to a true NLO correction to Z + jet production, can be separated by employing a veto against 
a second hard jet in real-emission events. Hence, we present NLO QCD corrections with a jet 

veto ( 5 QCD^eto> ^QCD.veto) and without a jet veto (^cd Z > <*qcd)- 

Using a jet veto based on a fixed pt value for the second jet is not well suited. It will either 
cut away relatively collinear emission events in the high-px tails of the leading-jet distribution 
(leading to large negative corrections) or it has to be chosen too large to be effective in the 
intermediate-pr parts of the distribution. Hence, building on our experience from W + jet 
production, we apply the jet veto defined at the end of Section 13.2.21 We have checked that this 
jet veto indeed effectively removes events with back-to-back jet kinematics. 

We also investigate the impact of the photon-induced tree- level processes (|2.5p and (|2.6p . 
Since the LO photon-induced cross section is a small effect, we show its relative impact <5 7 with 
respect to the LO cross section at 0(a 2 a s ) where initial states with photons are not taken into 
account. The NLO QCD corrections to these channels as well as the interference contributions 
discussed at the end of Section 12.31 are neglected in the following. 

In Tables [THSwe show the LO integrated cross sections, the corresponding relative corrections 
introduced above, and the NLO cross section Of^/^^ including the EW corrections for bare 
muons, the photon-induced processes, and the QCD corrections with the jet veto for the variable 
scale choice for different cuts on the transverse momentum of the leading jet and the dilepton 
invariant mass. All other cuts and the corresponding event selection follow our default choice as 
introduced in Section [3.2i In Figures [614131 we show for various observables the LO distribution 
and the distribution including the full set of corrections, i.e. EW corrections <5ew 5 the contribution 
of the photon-induced processes <5 7 , and the QCD corrections <5qcd- The various contributions 
to the corrections are also shown separately relative to the LO. All results are discussed in detail 
in the following subsections. 

3.3.1 Transverse momentum of the leading jet 

Tables HH3] show the LO predictions and the above corrections for different cuts on the of 
the leading jet prjet- All integrated cross sections and, hence, the corrections are dominated by 
events close to the lowest accepted pr,jet> as can be seen by the rapid decrease of the integrated 
cross section when increasing the PTjet cut. 

For the most inclusive cross sections (left columns in the tables) the EW corrections are 
at the percent level and negative. With increasing pr,jet) the relevant CM energies rise, and 
the well-known Sudakov logarithms in the virtual EW corrections start to dominate the total 
corrections as expected. For PTjet ~ 1000 GeV, the EW corrections are at the level of —25%. 
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pp -)■ l+l~ jet + X at = 14 TeV 



PTjet/ GeV 


25 - oo 


50-oo 


100 - oo 


200 - oo 


500 - oo 


1000 - oo 


^ =Mz /fb 


123491(7) 


44603(2) 


11364.4(5) 


1813.26(8) 


64.120(2) 


2.11859(6) 


aj a 7fb 


122024(7) 


43254(2) 


10445.0(4) 


1475.76(6) 


38.648(1) 


0.90847(3) 




-4.2 


-4.5 


-5.1 


-8.5 


-17.4(1) 


-27.0(1) 


rrec ,var /ny 
°EW / /o 


-2.8 


-3.2 


-4.2 


-7.8 


-16.7(1) 


-26.3(1) 




35.8(1) 


48.7(1) 


63.9(1) 


86.9(1) 


142.6(1) 


210.5(1) 




35.9(1) 


50.1(1) 


70.0(1) 


107.0(1) 


217.3(1) 


403.8(1) 


^QCD.veto/^ 


13.1(1) 


9.8(1) 


14.1(1) 


13.4(1) 


-2.4(1) 


-29.1(1) 


^QCD,vcto/% 


14.1(1) 


12.6(1) 


21.9(1) 


32.3(1) 


44.8(1) 


54.2(1) 




0.1 


0.2 


0.2 


0.4 


0.6 


1.0 


U + U~ ,var / ri 
<ullveto 


134266(49) 


46852(20) 


12223(4) 


1832.5(8) 


49.45(2) 


1.1649(7) 



Table 1: Integrated cross sections for different cuts on the pr of the leading jet (jet with 
highest pr) at the LHC with y/s = 14 TeV. We show the LO results both for a variable 
and for a constant scale. The relative EW corrections 5ew are given with and without 
lepton-photon recombination. The QCD corrections Sqcd are presented for a fixed as well 
as a for variable scale and with or without employing a veto on a second hard jet. The EW 
corrections and the corrections due to photon- induced processes, <5 7 , are presented for the 

variable scale. Finally, we show the full NLO cross section Of^i/^™^ The error from the 
Monte Carlo integration for the last digit (s) is given in parenthesis as far as significant. 
See text for details. 



This behaviour is generic and also present in all other observables where the cross section is 
dominated by events with high CM energies. 

It is evident from Tables [1] and [2] that the CM energy of the LHC plays a minor role for the 
EW corrections, in particular for less restrictive cuts. Since the transverse momentum of the 
leading jet depends only indirectly on the treatment of the lepton-photon system the corrections 
with and without recombination only slightly differ. The results for the corresponding differential 
distribution are shown in Figure [6] for the LHC with y/s = 14 TeV and the Tevatron. 

The transverse-momentum distribution of the jet has been calculated using the approximation 
of a stable, on-shell Z boson in Refs. [30, 31J. Of course, in contrast to our calculation, in on-shell 
Z+jet production it is not possible to apply various event-selection cuts to the leptonic final state, 
because the degrees of freedom related to the decaying Z boson are implicitly integrated out. 
Nevertheless, the relative EW corrections at high momentum transfer are dominated by Sudakov 
logarithms of the form In 2 (s/M|) that give rise to large process-independent contributions which 
factorize from the LO cross section. Therefore, they are expected to result in a similar behaviour 
for both the on- and off-shell corrections. Comparing our results for the leading-jet PTjet ha 
Figure [6] with Figure 5 in Ref. [31 j . we find agreement within 1—2% for PTjet > 200 GeV. Only 
at smaller prjetj the details of the event definition start to be more relevant and the on-shell 
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pp -> l+l- jet + X at = 7TeV 



PTjet/ GeV 


25- oo 


50-oo 


100 -oo 


200 - oo 


500-oo 


1000 - oo 


cj(;- Mz /fb 

aS ar /fb 


53029(3) 
51949(3) 


17736.0(5) 
16881.3(5) 


3939.51(9) 
3482.85(6) 


471.85(1) 
357.071(7) 


7.4538(2) 
3.90038(8) 


0.06464(6) 
0.02139(2) 


_l_ _ 

<W A 


-4.2 
-2.7 


-4.4 
-3.1 


-4.9 
-4.0 


-8.0 
-7.3 


-16.6(1) 
-15.9(1) 


-26.0(1) 
-25.2(1) 


<Sqcd z /% 


35.1(1) 
36.3(1) 


43.3(1) 
47.0(1) 


52.2(1) 
63.9(1) 


68.9(1) 
100.9(1) 


119.2(1) 
229.7(1) 


191.6(1) 
505.0(2) 


r/i=M z /C w 

°QCD,veto/ /0 
^QCD,veto/% 


15.4(1) 
17.2(1) 


8.7(1) 
13.4(1) 


8.9(1) 
20.7(1) 


4.0(1) 
30.9(1) 


-15.3(1) 
48.2(1) 


-44.9(1) 
69.4(1) 




0.2 


0.2 


0.3 


0.4 


0.8 


1.9 


u~^~ u~ ,var / n 
^fullvcto /fb 


58823(23) 


18438(8) 


4045(2) 


440.2(2) 


5.168(2) 


0.03106(2) 



Table 2: Integrated cross sections for different cuts on the px of the leading jet at the 
LHC with yfs = 7TeV. See caption of Table [Hand text for details. 



pp -»■ l+l~ jet + X at y/s = 1.96 TeV 



PTJet/ GeV 


25 - oo 


50-oo 


75 — oo 


100 -oo 


200 - oo 


300 - oo 


a$= Mz /fb 


9648.3(3) 


2440.3(1) 


840.01(5) 


340.51(2) 


17.679(1) 


1.3935(1) 


0-8*7 fb 


9365.0(4) 


2268.2(1) 


741.47(5) 


284.88(2) 


12.0551(9) 


0.79178(7) 




-4.1 


-4.2 


-4.0 


-4.2 


-6.3 


-8.4 




-2.6 


-2.8 


-2.9 


-3.3 


-5.5 


-7.7 


^cd7% 


30.7(1) 


25.7(1) 


20.3(1) 


14.2(1) 


-8.7(1) 


-31.0(1) 


w% 


32.9(1) 


31.9(1) 


31.9(1) 


31.8(1) 


33.4(1) 


33.8(1) 


^QCD.veto/^ 


19.8(1) 


5.8(1) 


0.4(1) 


-5.6(1) 


-29.0(1) 


-50.6(1) 


^QCD,veto/% 


22.5(1) 


12.4(1) 


12.3(1) 


11.4(1) 


9.1(1) 


7.1(1) 


6™/% 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


,vax , « 
°full,veto / 10 


11102(6) 


2458(2) 


804(1) 


306.2(3) 


12.43(1) 


0.784(1) 



Table 3: Integrated cross sections for different cuts on the pt of the leading jet at the 
Tevatron. See caption of Table [Hand text for details. 
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Figure 6: LO and fully corrected distribution (top), corresponding relative EW and 
photon-induced corrections (middle), and relative QCD corrections (bottom) for the trans- 
verse momentum of the leading jet at the LHC (left) and the Tevatron (right). See text 
for details. 
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results deviate as expected. However, the good agreement between the two calculations is still 
remarkable, having in mind that in Ref. [31] both the QED corrections have been completely 
neglected and a different renormalization scheme (MS) has been adopted. 

The contribution <L from the photon-induced processes are small and only reach up to a few 
percent for large cut values where the EW and QCD corrections to the dominating tree processes 
are by far larger. Hence, we can safely neglect the corresponding NLO QCD corrections which 
are formally of the same order as the EW corrections to the partonic processes without a photon 
in the initial state. 

The qualitative features of the corrections at the Tevatron are very similar to those at the 
LHC. Of course, at the Tevatron the high-energy (Sudakov) regime is not as accessible as at the 
LHC, but the onset of the Sudakov dominance is nevertheless visible as can be seen in Table [3]and 
Figure El We have adapted the range for the different integrated cross sections to the kinematic 
reach of the Tevatron. 

Turning to the NLO QCD results at the LHC, we observe exactly the same qualitative results 
found for W + jet production in Ref. [29J . As discussed above, the differential cross section for 
large PT,jet> as shown in Figure O contains large contributions from a completely different class 
of events for which two jets recoil against each other. Hence, the corrections are huge. The 
correction 5q^ z is smaller than 5qcd because it is defined relative to a larger LO cross section. 
In absolute size, however, the NLO corrections are similar. Using the jet veto proposed at the end 
of Section 13.2.21 the corrections are reduced, and 5qq D rises only to the 50% level for large cut 

values. The fixed scale choice accidentally leads to even smaller corrections 8q^ z . As expected, 
the large discrepancy of the LO results for the two scale choices is largely removed by including 
the NLO corrections. However, varying the exact definition of the jet veto, the variable scale 
turns out to be more robust. 

For the smaller CM energy ^/s = 7 TeV the same qualitative behaviour is found (see Table [2] 
for the integrated results with varying cuts). At the Tevatron, the jet veto is not as important 
as at the LHC because of its kinematical limitations. On the other hand, as expected, the fixed 
scale choice leads to more and more negative corrections with increasing prjet ; in particular when 
employing a jet veto. Using the variable scale stable results are obtained. 

3.3.2 Dilepton invariant mass 

Tables 0HS1 show the analogous results for a variation of cuts on the dilepton invariant mass 
Mi[. The corresponding differential distributions for the Tevatron and the LHC at 14 TeV are 
displayed in Figure [71 Again, all other cuts and the corresponding event selection follow our 
default choice as introduced in Section 13.21 

In Tables HHSl only the cross section in the first column includes resonant Z production. All 
the less inclusive cross sections only contain the tail of the distribution generated by (far) off-shell 
Z bosons and photons. Hence, the cross sections are smaller than for comparable cuts on prjet- 
Note that for a large cut on Mn the cross section is dominated by events with a CM energy 
close to the Mu cut since the leptons are mainly produced back to back and the additional jet is 
relatively soft. In comparison, for a large cut on PTjet the CM energy has to be more than twice 
the cut value because the transverse momentum of the jet has to be balanced. 

The dilepton invariant-mass distribution is a crucial observable at the LHC, in particular 
for detector calibration and also in searches for heavy dilepton resonances. For calibration, the 
Z+jet channel is the main source for boosted Z bosons close to their mass shell resulting in high- 
energy leptons which are balanced by the hard jet. Hence, the precise theoretical understanding 
of the Mu distribution, in particular around the Z-boson peak, is mandatory. 



16 



pp -> l+r jet + X at = 14TeV 



Mu/GeV 


50 — OO 


100 — oo 


200 — oo 


500 — oo 


1000 — oo 


2000 — oo 




123491(7) 


7696.9(8) 


628.47(6) 


49.380(6) 


5.1124(6) 


0.27096(3) 


cr /tb 


lzzUz4( f j 


tDoo.z{o ) 


DUz.4o(0 J 


40. /OU(0 J 


4.oyiy(b J 


U.zo4oo\o ) 


var,Mn / ™ 

(Tq /ID 


1z18o8( r ) 


1 A in Q iQ\ 

/419.8(8J 


r on *VA ( K\ 

539. /4(5J 


Q A i no ( A \ 
34.1Uz(4) 


z. /y58(4j 


0.1U831(1 j 


<f ,va 7% 


-4.2 


-9.3(1) 


-5.7 


-9.5 


-15.1(1) 


-23.8(1) 


t-rcc .var /fry 


-2.8 


-5.2 


-3.0 


-5.8 


-10.3(1) 


-17.1(1) 




35.8(1) 


28.9(1) 


12.0(1) 


-11.3(1) 


-34.4(1) 


-62.7(1) 




35.9(1) 


29.7(1) 


14.7(1) 


-5.7(1) 


-25.8(1) 


-50.8(1) 




36.1(1) 


30.8(1) 


24.7(1) 


23.6(1) 


25.9(3) 


31.4(3) 


°QCD,vcto/ /0 


13.1(1) 


D.o^l ) 


-9-5(1) 


— oz.y^i ) 


— OO.o^l ) 


— o0.o(^l J 


^QCD,veto/% 


14.1(1) 


8.7(1) 


-5.4(1) 


-25.4(1) 


-46.0(1) 


-71.3(1) 


cVa^M;; /C r/ 

°QCD,veto/ 70 


14.3(1) 


10.7(2) 


7.4(1) 


8.0(1) 


10.9(3) 


16.7(3) 




0.1 


0.9 


2.7 


2.9 


2.6 


2.3 


u+ M - ,var / r, 
"full.veto / ID 


134266(49) 


7580(9) 


551.9(4) 


31.10(5) 


1.906(4) 


0.0167(3) 



Table 4: Integrated cross sections for different cuts on the dilepton invariant mass M\\ 
at the LHC with yfs = 14TeV. In addition to the cross sections and corrections given in 
Tables HH3] we also show results for an alternative variable scale choice introduced at the 
end of Section 13.3.21 See text for details. 

Concerning the EW corrections, the large Sudakov logarithms again result in large negative 
corrections at large values of the Mu cut. They reach the level of —20% for Mn > 2000 GeV. 
However, they cannot be compared to earlier on-shell computations since only off-shell Z bosons 
and photons contribute in the high-M// region as mentioned above. 

As expected, the corrections for bare muons are larger since photons, being radiated 
collinearly to one of the charged leptons, carry away momentum and reduce the dilepton in- 
variant mass if not recombined with the emitting lepton. This effect can be observed for any cut 
value in Tables HHS1 However, it is most prominent around the Z-boson peak in the differential 
distribution displayed in Figure [7J For bare muons, events that are enhanced by muon-mass log- 
arithms are shifted to lower bins in the distributions. Hence, the correction is large and negative 
up to —20% around the Z peak. The correction below the Z resonance is positive. While the 
absolute correction compared to the differential cross section at the peak is relatively small, the 
relative correction to the distribution is huge because the LO result is relatively small away from 
the narrow Z-boson resonance. It amounts to more than 100% around a dilepton invariant mass 
of 75 GeV. The large correction, however, does not signal a breakdown of perturbation theory 
but only reflects the particular kinematic situation around the Z peak. If photons in a small cone 
around a radiating lepton are recombined, the corrections from almost collinear photons outside 
the cone are still sizeable and amount to roughly half the size of the corrections for bare muons. 
The size of the corrections is further reduced if the cone size for the lepton-photon recombina- 
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pp -> l+r jet + X at yfs = 7TeV 



M u /GeY 


50-oo 


100 -oo 


200 - oo 


500 -oo 


1000 - oo 


2000 - oo 


Uq z /fb 
^ ar /fb 

a var,A/ n/fb 


53029(3) 
51949(3) 
51865(3) 


3259.1(4) 
3171.4(4) 
3064.5(4) 


242.06(3) 
230.25(2) 
194.99(2) 


13.300(2) 
12.283(2) 
8.235(1) 


0.7276(1) 
0.65366(9) 
0.33869(5) 


0.009264(9) 
0.008059(8) 
0.002829(1) 


^-,var /% 


-4.2 
-2.7 


-9.3 
-5.2 


-5.8 
-3.0 


-10.2 
-6.1 


-16.7(1) 
-11.0(1) 


-27.7(1) 
-19.3(1) 


Cd z /% 

[»var A -/j j / f>y 

A v °' i i lvJ -LL j U7 

°QCD / /o 


35.1(1) 
36.3(1) 
36.5(1) 


28.3(1) 
30.4(1) 
32.6(1) 


11.0(1) 
14.8(1) 
29.4(1) 


-11.8(1) 
-5.3(1) 
32.4(2) 


-34.0(1) 
-25.0(1) 
38.9(3) 


-63.5(7) 
-50.4(2) 
53.0(3) 


A /i=M z /C / 
°QCD,veto/ /0 

^QCD,veto/% 
<-var,M n /0 / 
°QCD,veto/ /0 


10.4(1 J 

17.2(1) 
17.6(1) 


9.3(1) 
12.0(1) 
15.3(1) 


-0.8(1) 
-2.1(1) 
14.7(1) 


on r\(-\ \ 
— oU.U(l ) 

-22.2(1) 

19.1(2) 


— OZ. 

-42.3(1) 
26.1(3) 


— oi.oyl ) 

-67.1(1) 
41.0(3) 


^ ar /% 


0.2 


1.0 


2.9 


3.0 


2.8 


3.0 


u + u,~ ,var / n 
<llvcto /ft 


58823(23) 


3290(4) 


218.7(1) 


8.67(1) 


0.2861(6) 


0.00065(1) 



Table 5: Integrated cross sections for different cuts on the dilepton invariant mass M\\ 
at the LHC with sfs = 7TeV. See caption of Table H] and text for details. 



tion is increased. To reach an accuracy level of a few percent near the resonance, multi-photon 
radiation should be included in the calculation. 

The CM energy at the LHC hardly plays any role for the size of the relative EW corrections. 
In Figure [8] we show the LO cross section and EW corrections for y/s = 7TeV. Only the LO 
cross section strongly depends on the CM energy. Even the results for the EW corrections at 
the Tevatron (see Figure [7]) hardly differ from the ones at the LHC. 

Since the dilepton mass is a property of the lepton-photon system only and it is invariant 
under boosts when recoiling against a hard QCD jet, our results can be compared to the EW 
corrections to the line shape of inclusive Z-boson production as, for instance, recently investigated 
in Ref. [13] . It turns out that the EW corrections are quite similar around the Z-boson peak. 
However, in the lower tail of the distributions the corrections only reach up to 80% in the inclusive 
analysis (see e.g. Figure 12 in Ref. [13]). A large part of this deviation can be attributed to the 
fact that the LO result for Z + jet production drops off faster than in the inclusive case, so that 
the relative corrections are different. 

Turning again to NLO QCD, the corrections around the Z resonance in the invariant-mass 
distribution of the leptons are pretty flat since QCD can only indirectly affect the line shape of 
the colour-neutral Z boson. Depending on the exact event definition, the corrections amount to 
a few tens of percent. However, as shown in Tabled! i n t ne high-invariant-mass tail both scale 
choices fail to reflect the kinematical situation, since the production of a far off-shell Z boson is 
dominated by the region near the threshold set by the cut on Mu. In this region the Z boson 
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Figure 7: LO and fully corrected distribution (top), corresponding relative EW and 
photon-induced corrections (middle), and relative QCD corrections (bottom) for the dilep- 
ton mass at the LHC (left) and the Tevatron (right). 
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Table 6: Integrated cross sections for different cuts on the dilepton invariant mass Mu 
at the Tevatron. See caption of Table 2] and text for details. 




Figure 8: LO and fully corrected distribution (left) and corresponding relative EW and 
photon-induced corrections (right) for the dilepton mass at the LHC with y/s = 7 TeV. 
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decays mainly to back-to-back leptons with relatively soft jet activity. Hence, Sq^ z as well as 
Sqqd become large and negative. Consequently, a variable scale choice based on the invariant 
mass of the lepton pair 



var,M„ = w,M„ = ^var,M„ = I + (p had )2 (3.5) 



reflects the underlying kinematics in a better way. Thus, for this particular observable, we also 
present LO and NLO QCD predictions for this scale choice in Tables 0HH1 Indeed, a significant 
stabilization of the relative QCD corrections, especially for large Mu cuts at the LHC, is found. 
The relative EW corrections do only change insignificantly when using the alternative scale. 



3.3.3 Transverse momentum of the charged leptons 

In contrast to the invariant mass of the two charged leptons the transverse momentum of each of 
the leptons is sensitive to the recoil due to the hard jet in the event. Hence, the Jacobian peak 
of the transverse momentum at half the Z-boson mass in the inclusive Z-boson sample is washed 
out when an additional hard jet is present. The LO results and the corresponding correction for 
the positively charged lepton are shown in Figure [9l While the corrections are equivalent for the 
negatively charged lepton at the Tevatron, at the LHC they differ at the percent level but are 
qualitatively similar. 

The EW corrections are at the percent level and quickly increase in size in the tail of the 
distribution for p T l+ > 60 GeV, in particular for bare muons where they almost reach —10%. For 
even higher Pt,z+ > 100 GeV the typical increase of the correction due to the Sudakov logarithms 
can be observed. However, in contrast to the transverse-momentum distribution of the leading 
jet, the lepton-photon recombination still matters and the corrections differ by roughly 5% in 
the whole tail of the distribution. 

The QCD corrections 5qcd for Ptj+ show pronounced dips where the LO cross section has 
peaks (see Figured]). The real corrections do not particularly populate the regions of the dis- 
tributions that are enhanced due to the particular LO kinematics. For f>iy+ < 100 GeV, <Sqcd Z 
and Sqqy) are practically identical. However, for larger transverse momenta, the corrections for 

the two scale choices differ significantly. Here, Sq^ z grows large and negative to compensate 
for the overestimated LO cross section. This is expected, since the hard jet recoiling against the 
high-pT lepton should be reflected in the scale choice. 

At the Tevatron, the kinematical features are less pronounced and the deviation using a fixed 
scale starts already at pr,z+ ~ 70 GeV. However, all qualitative features for the EW as well as 
the QCD corrections are the same. 



3.3.4 Transverse mass of the charged leptons 

For Z-boson production, the transverse mass does not play a central role because the lepton 
system can be fully reconstructed. However, it is instructive to compare the transverse-mass 
distribution for Z+jet production with the corresponding distribution for W + jet. For W + jet, 
due to the neutrino in the final state, only the transverse mass but not the invariant mass of 
the two leptons can be measured. The distribution is displayed in Figure [TUl It shows all the 
features of the analysis for W + jet discussed in detail in Ref. |29j . However, the EW corrections 
are roughly a factor of two larger and reach the level of —20% for bare muons in the interesting 
region around the Jacobian peak, simply because there are two charged leptons in the final state 
which radiate photons. This example indicates that one has to be especially careful when using 
ratios of W and Z cross sections with the aim to reduce theoretical and parametric (such as 
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Figure 9: LO and fully corrected distribution (top), corresponding relative EW and 
photon-induced corrections (middle), and relative QCD corrections (bottom) for the trans- 
verse momentum of the positively charged lepton at the LHC (left) and the Tevatron 
(right). 



22 



PDF) uncertainties. Even though the two systems are almost equivalent from the point of view 
of QCD, EW corrections differ significantly. 

3.3.5 Rapidity distributions 

In Figure [Til we display the results on the rapidity distribution for the positively charged lepton. 
For the LHC, the EW and QCD corrections are rather flat and only slightly increase the differ- 
ential cross section in the forward and backward regions. The distribution and the corrections 
are symmetric with respect to yi+ = 0. At the Tevatron, this symmetry is only slightly disturbed 
at LO and NLO QCD. However, the EW corrections are asymmetric, large and negative for 
Di+ = —2.5, and almost zero for yi+ = 2.5. Of course, the distribution for the rapidity of the 
negatively charged lepton shows the analogous behaviour with reversed rapidities w.r.t. the case 
of the l+. 

Concerning the rapidity of the leading jet yj e t, the EW corrections are flat and do not disturb 
the LO shapes of the distribution, as can be seen in Figure H2J The QCD corrections are only 
slightly larger in the forward and backward regions compared to the central region at the LHC. 

Since both leptons can be fully reconstructed, also the rapidity yu of the final-state lepton 
pair — which resembles the rapidity of the intermediate boson — is experimentally accessible. The 
corresponding results are shown in Figure [I3j 

4 Conclusions 

Following our study on W+jet production [29], we have presented the first calculation of the full 
EW NLO corrections to the production of two opposite-sign charged leptons in association with 
a hard jet at hadron colliders. For many observables the cross section is dominated by on-shell 
Z + jet production with a subsequent leptonic Z-boson decay. However, in our calculation all 
off-shell effects as well as the contributions of and the interference with an intermediate photon 
are taken into account, so that also observables that are not dominated by on-shell Z bosons are 
described with NLO accuracy in view of both QCD and EW corrections. 

We have implemented our results in a flexible Monte Carlo code which can model the ex- 
perimental event definition at the NLO parton level. The distinction of Z -(-jet and Z + photon 
production is consistently implemented by making use of the measured quark-to-photon frag- 
mentation function. We have also recalculated the NLO QCD corrections supporting a phase- 
space-dependent scale choice. Photon-induced processes are included at leading order but turn 
out to be phenomenologically unimportant. 

The presented EW corrections are particularly large for the Z-boson line shape, i.e. the 
dilepton-invariant-mass distribution, mainly due to collinear final-state radiation. At the peak of 
the distribution, the corrections reach —20% for bare muons while they are at the order of 100% 
in the lower tail of the distribution, even larger than for the inclusive Z-boson line shape where 
no additional hard jet in the final state is demanded. When the EW corrections are not enhanced 
by final-state radiation due to particular kinematics, they are typically negative and at the level 
of a few percent. However, in the tail of distributions which are sensitive to large CM energies, 
the well-known Sudakov logarithms become dominant, and the EW corrections increase up to 
—25% at partonic \fs ~ 2 TeV. For the pt distribution of the jet, these results agree with earlier 
results in the on-shell approximation for the Z boson [31J. The QCD corrections have a typical 
size of a few tens of percent. However, they can become extremely large (hundreds of percent) at 
large jet pt unless a sensible veto on a second hard jet is applied. The presented integrated cross 
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Figure 10: LO and fully corrected distribution (top), corresponding relative EW and 
photon-induced corrections (middle), and relative QCD corrections (bottom) for the trans- 
verse mass of the two charged leptons at the LHC (left) and the Tevatron (right). 
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Figure 12: LO and fully corrected distribution (top), corresponding relative EW and 
photon-induced corrections (middle), and relative QCD corrections (bottom) for the ra- 
pidity of the leading jet at the LHC (left) and the Tevatron (right). 
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sections and differential distributions demonstrate the applicability and flexibility of our setup 
which will be useful for accurate predictions of any observable in the investigated final state. 

The importance of the neutral-current Drell-Yan process as a testing ground for perturba- 
tive calculations and for understanding and calibrating the detectors at the LHC can hardly be 
overestimated. Our calculation extends the availability of theoretical predictions for this pro- 
cess class to the EW corrections to associated production of Z bosons with a hard jet. Our 
calculation also constitutes one important part of a full NNLO prediction of the mixed EW and 
QCD corrections for inclusive Z production. In the range of intermediate and large transverse 
momenta of the additional hard jet our calculation delivers state-of-the-art predictions, for small 
transverse momenta the pure NLO calculation should of course be improved by dedicated QCD 
resummations, a task that goes beyond the scope of this paper. 
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